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Abstract 



Cosmological TV-body simulations are now being performed using Newtonian gravity on scales 

X> 

Q ' larger than the Hubble radius. It is well known that a uniformly expanding, homogeneous ball of 

■ dust in Newtonian gravity satisfies the same equations as arise in relativistic FLRW cosmology, 



and it also is known that a correspondence between Newtonian and relativistic dust cosmologies 
continues to hold in linearized perturbation theory in the marginally bound/spatially flat case. 
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Nevertheless, it is far from obvious that Newtonian gravity can provide a good global description of 
■ an inhomogeneous cosmology when there is significant nonlinear dynamical behavior at small scales. 

We investigate this issue in the light of a perturbative framework that we have recently developed 



Q. 



CTV which allows for such nonlinearity at small scales. We propose a relatively straightforward 

"dictionary" — which is exact at the linearized level — that maps Newtonian dust cosmologies into 
general relativistic dust cosmologies, and we use our "ordering scheme" to determine the degree to 
which the resulting metric and matter distribution solve Einstein's equation. We find that, within 
^ ■ our ordering scheme, Einstein's equation fails to hold at "order 1" at small scales and at "order e" at 

large scales. We then find the additional corrections to the metric and matter distribution needed 
to satisfy Einstein's equation to these orders. While these corrections are of some interest in their 
own right, our main purpose in calculating them is that their smallness should provide a criterion 
for the validity of the original "dictionary" (as well as simplified versions of this dictionary). We 
expect that, in realistic Newtonian cosmologies, these additional corrections will be very small; 
if so, this should provide strong justification for the use of Newtonian simulations to describe 
relativistic cosmologies, even on scales larger than the Hubble radius. 
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I. INTRODUCTION 



With the improvements in computational abilities that have taken place in recent years, 
it is now feasible to do numerical simulations of structure formation in cosmology on scales 
comparable to — or even larger than — the Hubble radius. Such simulations are being carried 
out by a number of groups However, these simulations are being carried out using 

Newtonian gravity. Although it would appear reasonable to expect Newtonian gravity to 
yield a good description of phenomena on scales much smaller than the Hubble radius — 
except, of course, in the immediate vicinity of strong field objects — at first thought, it 
might seem absurd that it could be expected to yield a reliable description of phenomena 
on scales comparable to, or larger than, the Hubble radius. After all, Newtonian gravity 
posits forces that act instantaneously over arbitrarily large distances, whereas the dynamical 
evolution laws of general relativity assert that all influences propagate causally and that the 
distribution of matter outside of one's past light cone is irrelevant. Similarly, the Newtonian 
gravity description of the Hubble expansion involves relative motion of bodies, whereas the 
general relativistic description involves the expansion of space. Why should Newtonian 
gravity give an accurate description of behavior on scales comparable to — or greater than — 
the Hubble radius, when the relative velocity of bodies is comparable to — or greater than— 
the speed of light? 

n 

Nevertheless, as we shall review in the next section, it is well known (see, e.g., |6[) that 
under the assumptions of spatial homogeneity and isotropy, the equations for a uniformly 
expanding pressureless fluid ("dust") in Newtonian gravity are identical to the dynamical 
equations for a dust filled Friedmann-Lemaitre- Robinson- Walker (FLRW) universe in general 
relativity — even in the case of nonvanishing spatial curvature. An explanation for this 
remarkable correspondence can be found from the fact that in both Newtonian gravity 
and general relativity, in the presence of spherical symmetry, the behavior of a co-moving 
ball of dust does not depend upon the distribution of matter outside of the bal{j} Thus, 
in both Newtonian gravity and general relativity, the dynamical behavior of a co-moving 
ball of dust in a homogeneous, isotropic universe is the same as it would be if that ball 
were placed in an empty, asymptotically flat spacetime. However, for a sufficiently small 
ball of dust in an otherwise empty spacetime, Newtonian gravity should be an excellent 
1 This fact is closely related to the fact that there is no gravitational field / curvature inside a spherical shell 
of matter in Newtonian gravity (by Newton's theorem) and general relativity (by Birkhoff's theorem). 
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approximation to general relativity. Thus, for a sufficiently small ball, the density and co- 
moving radius of the ball must satisfy the same dynamical equations in Newtonian gravity 
and general relativity. This implies that the equations for the density and comoving radius 
in a homogeneous, isotropic Newtonian dust cosmology must coincide with the equations for 
the density and scale factor in a FLRW dust cosmology. This correspondence continues to 
hold in the presence of a cosmological constant term in the Newtonian and general relativistic 
equations. 

The above argument relies crucially on exact spherical symmetry. Thus, one might ex- 
pect that no such correspondence between Newtonian and relativistic cosmologies would 
hold if one perturbs the homogeneous, isotropic solutions away from spherical symmetry. 
Remarkably, however, the correspondence between Newtonian and relativistic cosmologies 
extends into the regime of linearized perturbation theory in the case of perturbations off of a 
spatially flat FLRW dust cosmology. More precisely, as pointed out by Bardeen [7| and will 
be reviewed in the next section, the scalar gauge-invariant variables of linearized relativistic 
perturbation theory obey exactly the same equations as the variables describing linearized 
irrotational dust perturbations of the corresponding Newtonian cosmology. Furthermore, 
it is not difficult to see that this correspondence extends to the vector well, i.e., 

vector gauge-invariant variables of linearized relativistic perturbation theory obey exactly 
the same equations as the corresponding Newtonian variables describing vorticity perturba- 
tions. Thus, the scalar and vector sector^ of linearized relativistic perturbation theory off 
of a spatially flat FLRW dust model are in exact correspondence with arbitrary Newtonian 
perturbations off of the corresponding Newtonian dust cosmology. 

Further justification for the validity of the Newtonian approximation in cosmology is 
provided by the work of Oliynyk js, 9] (see also Futamase 10]). Oliynyk rigorously proved 
that for a given 3-torus Newtonian cosmology, there exists a one-parameter family of general 
relativistic solutions that limits to this Newtonian cosmology, thus showing that there are 
general relativistic solutions that are arbitrarily close to the Newtonian solution. However, 
for Oliynyk's one-parameter families, the ratio of the size of the universe to the Hubble radius 
goes to zero in the limitjfl. Thus, the general relativistic solutions proven by Oliynyk to be 

2 The tensor modes correspond to additional degrees of freedom present only in general relativity, and they 

have no Newtonian correspondence. 

3 Oliynyk formulated his limit as one in which the 3-torus remains of fixed size, but the speed of light — and 

the Hubble radius — goes to infinity. By re-scaling the spatial coordinates, his limit can be reformulated 
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very close to a Newtonian solution have size small compared with the Hubble radius, and thus 
have no "long wavelength part" . Thus, Oliynyk's results do not directly address the issue of 
whether Newtonian simulations on scales comparable to the Hubble radius correspond closely 
to a general relativistic solution, but it can be viewed as providing additional justification 
for the validity of Newtonian gravity on scales small compared with the Hubble radius. 

Taken together, the above considerations strongly suggest that for a universe that is suffi- 
ciently close to a spatially flat FLRW dust model, Newtonian gravity should provide a good 
description of structure formation on all scales. However, the situation is far from straight- 
forward for the following reasons: (i) Although, as described above, there is a correspondence 
at linearized order between Newtonian theory and general relativity, the "dictionary" needed 
to translate a linearized Newtonian solution into metric and matter perturbations in any 
particular gauge is nontrivial, and it is not obvious how this dictionary compares with stan- 
dard dictionaries used for the Newtonian and post-Newtonian approximations to general 
relativity on small scales. Thus, it is not obvious how to produce a "global dictionary" 
that works on all scales, (ii) If one has a candidate global dictionary, it is not obvious 
how to formulate criteria to determine whether the resulting general relativistic spacetime 
is "sufficiently close" to a solution to Einstein's equation to trust its predictions. The main 
complication here is that the failure to take post-Newtonian corrections into account on 
small scales will cause the general relativistic spacetime to fail to satisfy Einstein's equation 
by a larger amount than the failure to properly account in any way for the long wavelength 
perturbations. For most applications in cosmology, the tiny post-Newtonian corrections to 
the metric and matter motion on small scales are of no interest, but the leading order devia- 
tion of the metric and matter density from a FLRW model on large scales is of great interest. 
Thus, the proper criteria for being "sufficiently close" to a solution to Einstein's equation 
must take into account the distinction between small scales and large scales, (iii) One would 
like to know explicitly what the dominant corrections to Newtonian cosmology are, both to 
be able to quantitatively judge its reliability and to be able to make its predictions more 
accurate. 

The difficulties in addressing the above issues stem from the fact that the approximations 
of Newtonian gravity (which, a priori, is expected to be good on small scales) and linearized 
as one in which the speed of light remains constant and the Hubble radius goes to a well defined limit, 
but the size of the 3-torus then approaches zero in the limit. 



4 



perturbation theory (which, a priori, is expected to be good on large scales) are incompatible. 
Specifically, in the Newtonian gravity approximation certain nonlinear terms in the equations 
are kept (as they must be at small scales), but it is essential that time derivatives of quantities 
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12( | . By contrast, linearized perturbation 



be small compared with space derivatives js, 9, 
theory allows time derivatives of quantities to be comparable to their space derivatives (as 
they must be at large scales), but it is essential that all nonlinear terms be negligible. In 
order to properly treat phenomena on all scales, one needs an approximation scheme that 
can accommodate nonlinear phenomena on small scales but treats time derivatives on the 
same footing as space derivatives on large scales. We recently proposed an approach that 
accomplishes this [l|, and we will apply this approach here^. 

The main questions we wish to address in this paper can now be stated concretely as 
follows. Suppose that a Newtonian cosmological simulation has been performed on a 3-torus 
(i.e., periodic boundary conditions), where the size of the 3-torus may be larger than the 
Hubble radius. For convenience, we assume that the Newtonian solution has been presented 
as a continuum solution — i.e., that suitable smoothing has been done if the solution was pro- 
duced from an iV-body simulation. We would like to know the following: (1) What general 
relativistic spacetime and dust matter distribution should we associate to this Newtonian 
cosmology, i.e., what "global dictionary" should we use? (2) To what extent is this space- 
time a solution to Einstein's equation, i.e., what are the leading order terms in Einstein's 
equation that fail to be satisfied? (3) What are the leading order corrections to the metric 
and dust distribution that improve the accuracy of this solution, and how large are these 
corrections? 

Our approach will be to use the framework of |l| to provide a "counting scheme" for the 
sizes of terms in Einstein's equation. We will start with a candidate "global dictionary," 
which is suggested by the known correspondences between Newtonian gravity and general 
relativity in the exactly homogeneous and isotropic case and at the linearized level. We will 
then see that in our counting scheme, the resulting general relativistic spacetime fails to sat- 
isfy Einstein's equation to 0(1) at small scales and to 0(e) at large scales. The main effort 
in our paper will then be to find the corrections to the metric and dust distribution that, 
within our counting scheme, improve the accuracy of the solution to 0(1) at small scales 
and to 0(e) at large scales. It should be emphasized that we shall not prove existence of 



Our approach is closely related to 



13|; see also 



14]. 
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a one-parameter family of solutions to Einstein's equation with the properties we desire — a 
far more difficult task than solving for leading order corrections. Nevertheless, if the leading 
order corrections we obtain are small compared with terms appearing in the original global 
dictionary, we believe that this provides a strong indication that there is a general rela- 
tivistic solution that corresponds closely to the Newtonian cosmology. Conversely, if these 
corrections are not negligibly small, then either the Newtonian cosmology is not providing 
a sufficiently accurate representation of the general relativistic spacetime or the dictionary 
being used will have to be significantly modified. 

Our analysis also addresses concerns that have been expressed with regard to the use of 
a "Newtonianly perturbed FLRW metric," which corresponds to the using the "abridged 
dictionary" given by (I2.46p - fl2.48p below. Ishibashi and Wald [15] have argued that this 
metric should provide an excellent description of our universe. However, several authors 

nn 

[161 . 1171 ] have objected to the use of this metric on the grounds that, if taken literally, and 
dust peculiar velocities are ignored, then strong constraints relating to exact solutions of 
Einstein's equation apply, and the metric is only able to describe a spatially homogeneous 



continuum. Other concerns have been raised by Rasanen [18j. The spacetime metric and 
dust matter distribution that we produce in this paper — as summarized in section HV1 — solves 
Einstein's equation to a much higher degree of accuracy than the Newtonianly perturbed 
FLRW metric does, and, in particular, fully takes into account peculiar velocities and leading 
nonlinear terms in the Einstein equation. No inconsistencies of any kind are encountered 
in obtaining this much more accurate solution. Thus, the approximate solution considered 
in [fjj] should be fully justified provided only that the corrections to f l2.46p -f l2.48p given in 
section [IV] are negligibly small, as we argue is the case. 

We remark that if one has an equation S(F) = 0, one must draw a clear distinction 
between having a quantity / that approximately solves this equation (i.e., £(/) ~ 0) as 
compared with having a quantity / that is an approximate solution (i.e., / F for some 
exact solution S(F) = 0). If the equation is suitably well posed, if £(/) ~ 0, and if F is the 
exact solution with the same initial data as /, then / and F will remain close to each other for 
sufficiently early times. However, / may fail to remain close to F at late times because of the 
build-up of secular effects. For example, the Newtonian solution for the motion of Mercury 
solves the general relativistic equations of motion to an excellent approximation at all times, 
but provides a very poor approximation to the general relativistic solution for the position 



of Mercury after ~ 10 6 years. In this paper, we are concerned with the issue of obtaining 
general relativistic spacetimes that solve Einstein's equation to an excellent approximation 
at all times, but we will not be concerned with the issue of whether these spacetimes provide 
good global-in-time approximations to exact solutions of Einstein's equation. 

In the next section, we shall review the correspondence between homogeneous, isotropic 
Newtonian cosmology with dust matter and FLRW models in general relativity, as well 
as the correspondence at the linearized level between these models in the marginally 
bound/spatially flat case. On the basis of this correspondence, we will propose a dictionary 
fl2.40p -f l2.44p to translate Newtonian cosmologies into general relativistic spacetimes. In 
section IHH we will apply our counting scheme [l] to analyze how well Einstein's equation 
is being satisfied, and we will obtain corrections needed to satisfy Einstein's equation to 
0(1). We will then obtain the further modifications to the metric and dust distribution 
needed to obtain a solution to Einstein's equation to 0(e) at large scales. These corrections 
are of some interest in their own right. For example, as we shall see in Appendix [B], there 
are small modifications of some global properties of the cosmology, such as a slight modifi- 
cation of the expansion rate and the introduction an (even smaller) anisotropic expansion. 
However, our main purpose in determining these corrections is to provide a criterion for the 
validity of the Newtonian cosmology as translated into a general relativistic spacetime via 
the dictionary fl2T40D - fl2T44D and/or its abridgment f[2T46|) f[2T48|) or its simplification f|2T49]) 
(12.5ip : The full set of metric and matter corrections to our original dictionary are given 
by eqs. fl4.ip - fl4.13p and these corrections can be computed straightforwardly for any given 
Newtonian cosmology. The smallness of these corrections should provide a reliable criterion 
for judging the validity of using a Newtonian simulation with the dictionary fl2.40p - fl2.44p 
(or its abridgment or simplification) to describe a relativistic cosmology. 



II. BACKGROUND AND LINEARIZED CORRESPONDENCE 

In this section we shall review the correspondence between homogeneous, isotropic New- 
tonian dust cosmology and FLRW models, as well as the correspondence between linearized 
perturbations of these models. We shall then propose a dictionary — valid to linearized 
order — the accuracy of which will be evaluated and improved upon in the following section. 

In Newtonian gravity, the gravitational field is described by a Newtonian potential </>, 
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and the dust matter is described by a mass density p and a velocity field v % . The Newtonian 
field potential is related to the mass density by the Poisson equation 

<9'A0 + A = 47rp, (2.1) 

which we have generalized to allow for the presence of a cosmological constant A. In addition, 
the matter variables must satisfy mass conservation and Euler equations, which, for dust 
matter, take the form 

d t p + d i (pv i ) = 0, (2.2) 
dtiftf) + djipvV) = -p<9V . (2.3) 

In these equations, the flat Euclidean metric of space is used to contract indices. 



A. Background Correspondence 

As a cosmological ansatz, we seek a solution to the above equations of Newtonian gravity 
in which the density is spatially uniform, p = po(t), and the velocity field is uniformly 
expanding, v % = H{t)x % . Note that H is related to the radius, a, of any comoving ball by 

# = ~. (2.4) 

Since dtfj* = 3H(t), implies that 

d tPo + 3Hp = 0, (2.5) 

from which it follows that 

Po = Po,mitGT 3 . (2.6) 
Using mass conservation, we can eliminate p from the Euler equation, 

drf + v j d jV l = -<9> . (2.7) 

The Poisson equation, (12.1 p . has the non-singular solution 



2tt 2 A 
0o = ~^-por - - 

6 O 

Substituting for v and 0, we obtain 



Por 2 - -r 2 + A(t) . (2.8) 



dH 9 47r A . 
_ + ^ = __ P0 + _, (2.9) 



which is one of the Friedmann equations. To obtain the other Friedmann equation, we 
rewrite this equation as 

1 d 2 a 
a dt 2 

Integrating once, we obtain 

H 2 -- 

where k is & constant of integration. By choosing the size of the comoving ball appropriately, 
we may choose k to take the values 0, ±1. When A = 0, the value of k determines whether the 
universe is unbound and expands forever (k = — 1), is marginally bound and expands forever 
but with expansion velocity approaching zero (k = 0), or is bound and will recollapse within 
finite time (k = +1). Of course, in Newtonian gravity, k does not have any interpretation 
in terms of spatial curvature; space is always Euclidean. 

Equations (12. 9 h and ( 12. lip are precisely the equations satisfied by dust FLRW models in 
general relativity. The underlying reason for this exact correspondence was discussed in the 
Introduction. 



47T p ,i nit A 
3 a 3 3 



8vr A k 
~TPo + T ~ 3 > 



(2.10) 



(2.11) 



B. Linearized Correspondence 

We first re-write the exact Newtonian equations relative to some (arbitrarily chosen) 
background solution of the previous section. We introduce comoving coordinates 

x H = - . (2.12) 
a 

We then define the velocity variable v H by 

. • dx dx X da • j T _ ■ , . 

v n = a— = - -x l = v l - Ex 1 , 2.13 

dt' dt adt K ' 

so v n measures the velocity relative to the Hubble flow of the background solution. We also 
define density and potential deviations from the background, 5 and if), by 

p = Pa(l + 5), (2.14) 
= 0o + V- (2-15) 



In terms of these quantities, the Newtonian equations are 

&'di,i/) = 4:ira 2 po5, (2.16) 
d t ,8 + ~d v ((1 + <5K') = , (2.17) 

<W*' + -v'i'dfv* + Hv H ' = --&'ip . (2.18) 

From now on, since we will always work in comoving coordinates, we shall drop the primes. 
Next, we re-write these equations using "conformal time" r defined by 

dr 1 . 

— =-. (2.19) 

at a 

We also denote derivatives with respect to conformal time with overdots. In terms of the 

conformal time variable, the Newtonian background equations are 

a 2 8tc 2 A 
— = -z-a po + - 
a 2 3 3 



2 - o a z p + -a 2 -k, (2.20) 



d ( a\ An 2 A 2 

T" ~ =-^r«Po + T7« 2 , (2-21) 
dr \a J 3 3 

Po + 3-p = 0, (2.22) 
a 

and the Newtonian equations for the quantities describing the deviations from the back- 
ground are 

d%ip N = Ana 2 p 5 N , (2.23) 

6 N + d i ((l + 6 N )vi f )=0, (2.24) 

v* N + v j N d jV jf + -< = , (2.25) 

ct 

where we have now added a subscript A" so that these Newtonian quantities can be easily 
distinguished from the corresponding general relativistic quantities that we will introduce 
later. We emphasize that (I2.23p — (I2.25p are exact. We shall assume below that these equa- 
tions are solved on a 3-torus, i.e., a "box" at fixed comoving coordinates (of the background 
solution) with periodic boundary conditions. 

Linearizing (I2.23p -f l2.25p about the background solution, we obtain 





= 47ra 2 po4^ > 


(2.26) 


8$ + d lV ^ 


= 0, 


(2.27) 


. d (x)i 

V N + ~ V N 

(J/ 


= . 


(2.28) 
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We now compare the linearized Newtonian equations with the linearized general relativistic 
equations about a dust FLRW model. In [7], Bardeen decomposed linearized metric and 
stress-energy perturbations into their scalar, vector, and tensor parts, which evolve inde- 
pendently. He then introduced gauge invariant quantities to describe these perturbations. 
In the case of a perfect fluid, the two scalar gauge invariant variables describing metric per- 
turbations are related — in his notation, = — $a — so the scalar perturbations are fully 
described by (or the scalar part of the velocity perturbation, v l s , and a density 
perturbation variable e m (defined by eqs. (3.9)-(3.11) and (3.13) of [7J). We can similarly 
decompose a linearized Newtonian perturbation: ifjff and 5$ are scalar quantities, and the 
velocity perturbation can be decomposed as 

v { i* = «g? + t#» , (2.29) 

where can be written as a gradient and diV^ = 0. Newtonian perturbations have no 
tensor part. It is then straightforward to see that, as pointed out by Bardeen 3], under the 
correspondence 

</>£ } <— ► $4 = -$h , (2.30) 
«£? <— ► < , (2-31) 
6$ <— >> e m , (2.32) 

the linearized Newtonian equations become identical to the equations describing scalar per- 
turbations of a spatially flat dust cosmology as given by eqs. (4.3), (4.5) and (4.8) of 
Note that this correspondence holds only for perturbations of spatially flat models, i.e., there 
are additional terms in the linearized Einstein equation when the background solution has 
nonvanishing spatial curvature. 

It is not difficult to see that the exact correspondence between linearized Newtonian 
gravity and general relativistic perturbations of spatially flat models extends to vector per- 
turbations as well with 



(i)i 

' J Nv 



(2.33) 



where the gauge invariant quantity v l c is defined by Bardeen's eq. (3.23). Specifically, the 
vector part of the linearized Euler equation (I2.28j) for is identical to Bardeen's eq. (4.13) 
for i>*. There are no additional Newtonian equations for vector perturbations. However, there 
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is an additional general relativistic equation [Bardeen's eq. (4.12)] 



8>djtt x 




(2.34) 



which is a Poisson equation for a quantity \& l not present in Newtonian theory. On a torus, 
this equation has a solution if and only if there is no spatially homogeneous part of v % c , i.e., 



This equation must hold for all for all times r, but it is easily checked that this equation 
is preserved under time evolution, so it suffices to impose it at any one time. Thus, f l2.35|) 
is a constraint that must be imposed upon a linearized Newtonian solution in order that it 
correspond to a linearized solution of Einstein's equation under the correspondence f!2.33j) 

In summary, provided only that the constraint f |2.35|) is satisfied, there is an exact corre- 
spondence, given by f l2.30[) — ( ]2.34j) . between the complete linearized Newtonian equations for 
dust matter off of a homogeneous and isotropic background and the scalar and vector parts 
of the linearized Einstein equation off of a spatially flat dust FLRW background. As pre- 
viously noted, there are no counterparts to tensor perturbations in Newtonian gravity, i.e., 
general relativity has these additional degrees of freedom not present in Newtonian gravity. 

C. A Proposed Dictionary 

Based upon the results of the previous subsections, we now shall propose a dictionary 
that translates a solution (iPn,6n,v 1 n ) of the exact Newtonian equations (I2.23p - fl2.25p into 
a general relativistic spacetime metric g a b and dust matter stress-energy tensor T a b = pu a Ub- 
In the next section, we shall investigate the extent to which (g a b, T a b) satisfies Einstein's 
equation as well as what further corrections need to be made to (g a b, T ab ) to make it solve 
Einstein's equation to higher accuracy. 

First, the Newtonian equations fl2.23p - fl2.25p were written relative to a "background so- 
lution" of f)2.9p and f)2.1ip with "scale factor" a and mass density po- Since these Newtonian 
background equations are identical to the equations for a dust filled FLRW general rela- 
tivistic spacetime, we define our dictionary so that it associates the corresponding FLRW 
spacetime to this background solution. Thus, we have defined our dictionary for the case 
ipN = <$n = v l N = 0. In the following, we shall restrict consideration to the case where k = 




(2.35) 
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for the background solution, since this is the only case where we expect a good dictionary to 
exist when deviations from homogeneity and isotropy occur. We shall assume that periodic 
boundary conditions have been imposed on the Newtonian background solution, so that the 
corresponding FLRW background solution has 3-torus spatial slices. For convenience, we 
assume that the co-moving spatial coordinates of the Newtonian and FLRW background 
solutions range between and 1. 

In order for our dictionary to produce a definite general relativistic spacetime, we must 
make a choice of gauge for the metric. In the context of linearized perturbation theory, a 
natural and very useful gauge choice is the longitudinal gauge, in which the metric takes the 
form 

ds 2 = a 2 (r) [-(1 + 2A)dr 2 - 2B i dx i dr + ((1 + 2H L )5 ij + h i:j ) dx l dx j ] , (2.36) 

where d l Bi = 0, d^hij = and h % { = 0, and spatial indices i, j, k, . . . are raised and lowered 
with the background flat Euclidean metric 5ij. In the context of linearized perturbation 
theory, the quantities A, B i: H L , and hij represent the metric perturbation, and it can be 
shown that an arbitrary metric perturbation can be put in the form (12.36}) by an infinitesimal 
gauge transformation. It also can be shown that this gauge is essentially unique, i.e., there is 
essentially no additional gauge freedom that maintains the form (12. 36ft . However, linearized 
perturbation theory is not adequate for our purposes, since our dictionary is required to 
map Newtonian solutions that differ by a finite amount from the Newtonian background 
solution into metrics that differ by a finite amount from an FLRW model. Nevertheless, 
it should be possible to show via the implicit function theorem that for metrics that differ 
from an FLRW model by a sufficiently small but finite amount, the metric form (I2.36f) — 
with d l Bi = 0, d^hij = and h l i = — always can be imposed by a (nonlinear) gauge 
transformation. We shall not attempt to prove such a result here, and will merely take 
(I2.36P as an ansatz for the metric in constructing our dictionary. However, we believe that 
imposition of the metric form (I2.36|) does not involve any loss of generality if the metric is 
sufficiently close to an FLRW model. 

The stress-energy tensor of dust in the general relativistic spacetime takes the form 

Tab = pu a u b . (2.37) 

We define the three- velocity, v\ of the dust to be such that the components, it M , of the 
four-velocity in our gauge are proportional to (l,v l ). Normalizing using the metric form 
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(12.36j) . we obtain 



1 



= (!,««). 



(2.38) 



ay/1 + 2A + 2B jV J - ((1 + 2H L )6 jk + h jk ) v^v k 

Thus, this equation gives the formula for the 4-velocity u a appearing in (|2.37[) in terms of 
the 3-velocity v l that will be specified by our dictionary below. We define the fractional 
density perturbation 5 in the general relativistic model via 



p = Po{l + 6). 



(2.39) 



As already stated above, the Newtonian solution is specified by (iPn->$n, v n)- With the 
above gauge choice, the general relativistic spacetime and matter distribution is specified by 
(A, Bi, Hl, h^, S,v l ). Our proposed dictionary will therefore be defined by providing formulas 
for A, Bi, Hl, h^, 5, and v i in terms of the Newtonian variables. To obtain this dictionary, 
we start by taking the formulas that hold at linearized order under the correspondence of 
the previous section, which we obtain by expressing the Bardeen variables appearing in 
(12.301) — ( 0.331) in terms of (A,Bi,HL,hij,5,v l ). Then we improve our definition for Bi (and 
correspondingly for v l ) by requiring consistency with the nonlinear momentum constraint 
at small scales, leading to the replacement of p Q v l by p (l + 5n)v 1 . We thereby propose the 
following dictionary: 



A = -R L = i, 



N 



(i + s N y = (i + 5 N )(v i N + b 1 ) - (i + 6 N y N 

2 



5 = 5m - 



hij = , 



Anp a 2 



a \ . a ■ 

- ) Wn + ~Wn 

a J a 



and B % is the solution to the equation 



d 3 djB l 



wit 



-lQnp a 2 (Jl + c^K 
ljf| Bi = 0, where the overbar denotes spatial average, i.e., 



1 + S N )v l N 



(2.40) 
(2.41) 

(2.42) 
(2.43) 

(2.44) 



/= / a*xf. 



(2.45) 



The condition Bi = can be imposed by using the coordinate freedom x l 



•F*(t). 
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(Recall that the comoving spatial coordinates are assumed to range from to 1.) In ( 12.41 j) 
and (I2.44p . the notation | v denotes the "vector part" of a quantity in a decomposition of the 
type (J22U). 

At large scales, one would expect the vector part of v % to be negligible because in linear 
perturbation theory, vector modes are known to decay [7|. In addition, comparing the 
Poisson equation for B % with the Poisson equation for ip^, one would expect B l to be 
smaller than ip^ by order v/c at small scales, and thus B l should be negligible compared 
with A and Hl at all scales. Thus, (I2.4ip should yield a negligibly small correction to the 



equation v l 



'N- 



Thus, under normal circumstances, it should be acceptable to replace our 



proposed global dictionary with the following abridged version of the dictionary: 

A = -H L = i/j n , 



N > 



5 = 5 N - 



a 



(2.46) 
(2.47) 

(2.48) 



4irp a 2 

together with B % = = 0. This abridged dictionary corresponds to a continuum version 
of the dictionary^ given by Chisari and Zaldarriaga 19 ] . 

Finally, on small scales 5n should dominate the other terms appearing on the right side 
of (I2.42p . Thus, on scales much smaller than the Hubble radius, it should be possible to use 
the following simplified version of the dictionary: 



.4 



5 



5n , 



(2.49) 
(2.50) 
(2.51) 



together with B % = hij = 0. This simplified dictionary is very commonly assumed. However, 
on scales comparable to the Hubble radius, all terms in f|2.42p should be of comparable size, 
so if one is interested in investigating behavior on large scales, the full dictionary or abridged 
dictionary should be used. 

As explained above, our dictionary f l2.40p -f T2.44p has been constructed so as to produce 
a solution of Einstein's equation to linearized order in (A,B i ,HL,hij,5,v' t ). A Newtonian 
cosmology that corresponds to our universe should have ipN *C 1 and \v l N \ 1, but will 



6 Their definition of p differs from ours by a term involving the perturbed volume element. 
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normally have 5n S> 1 on small scales. Therefore, it is not obvious, a priori how accurate 
our dictionary is in producing a solution to Einstein's equation. In fact, it is clear that there 
may be difficulties in this regard because the dictionary should produce a spacetime that 
nearly satisfies the linearized Einstein equation on small scales, but the linearized Einstein 
equation is incompatible (via the linearized Bianchi identity) with the nonlinear dynamical 
behavior of matteiQ that occurs on small scales. We now investigate how close fl2.40p - fl2.44p 
comes to producing a solution to Einstein's equation. 

III. COUNTING SCHEME AND IMPROVED DICTIONARY 

As stated at the end of the previous section, we wish to determine how close our dictio- 
nary comes to producing a solution of Einstein's equation when the Newtonian cosmology 
has ipN "C 1 and \v l N \ <C 1, but may have 5n S> 1 at small scales. To analyze this is- 
sue, we need a consistent approximation scheme that can take advantage of the fact that 
the deviation of the metric from an FLRW model is small on all scales, but permits very 
large deviations of the stress-energy tensor from an FLRW model to occur on small scales. 
Such an approximation scheme was recently developed by us in [1] and used to analyze the 
backreaction effects of small scale inhomogeneities on large scale dynamics. We refer the 
reader to that reference for the precise mathematical formulation of the approximation. For 
our purposes here, it suffices to observe that, as in ordinary perturbation theory, in our 
approximation scheme there is a "small parameter" e (denoted A in [I]) that measures the 
deviation, = g a b — Qab \ °f the metric g a b from a background metric g^ b \ so ^ab = 0(e). 
However, unlike ordinary perturbation theory, first spacetime derivatives of 7 a b are allowed 
to be 0(1), and second spacetime derivatives of — and, hence, the deviations of the stress- 
energy tensor from the background stress-energy — are allowed to be 0(l/e). In particular, 
the quadratic products V dab^ fide and lab^ fide that appear in Einstein's equation are 
0(1), so our approximation scheme allows small scale inhomogeneities to affect the dynamics 
of the background metric. One of the main results of [l| is that, in fact, the only possible 
effect that these nonlinear terms can have on the dynamics of the background metric is to 
contribute an effective stress-energy that is traceless and has positive energy, corresponding 



7 In many references (see, e.g., 



19 



20j ). the linearized Einstein equation is written down together with 
the nonlinear dynamical equations for matter. This combined system of equations is mathematically 
inconsistent. 
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to the presence of gravitational radiation. For the present work, we assume that the universe 

contains a negligible amount of gravitational radiation, so that this effective stress-energy 

tensor can be set to zero, and the background metric (which has FLRW symmetry) therefore 

obeys the ordinary Einstein equation with dust stress energy tensor. 

In addition to analyzing the effects of small scale inhomogeneities on the dynamics of 

the background metric, in |l| perturbation theory was generalized to allow for significant 

nonlinearity at small scales, while at the same time maintaining a linearized description at 

large scales (see also [3|). In order to ascribe different behavior to perturbations at different 

scales these notions must of course be defined. In the notion of the "long wavelength 

part" of quantities was defined in a mathematically precise manner by considering the weak 

I I 

limit of these quantities as e — > 0. As explained in a ^ sufficiently small but finite e, 
this should correspond closely to taking an average over a spatial scalej L that is small 
compared with the background curvature (i.e., the Hubble radius) but sufficiently large that 
at this scale and beyond we have \S\ <C 1. For the present work, we shall identify the long 
wavelength part, A^:.. an , of a tensor field, A ai ... an , with the spatial average^ of its components 

4r-M» = (4*™#JC»0 = j d 3 x'W L (x-x')A^ n (x'), (3.1) 

using a suitable "window function," Wl{x — x'), of size L, i.e., a smooth function which is 
equal to 1 for a 2 \x — x'\ 2 < L, and which smoothly falls to outside of this regio ifl The 
requirement that L be much smaller than the background curvature scale ensures that this 
averaging process is well-defined, whereas the requirement that L be sufficiently large that 
\8\ <C 1 should ensure that the long wavelength parts of perturbations behave linearly. 
We define the "short wavelength part" of A ai ... an by 

thereby providing a decomposition of any quantity into its long and short wavelength parts. 
The framework of [l| allows one to make different assumptions in a mathematically consistent 
manner about the long and short wavelength parts of the various quantities. In particular, 



8 ForTlhe present universe, L ps 100 Mpc should meet these criteria. 

9 In 1] averages over (short) time intervals were also performed. The purely spatial averaging we perform 



here with a suitable "window function" corresponds precisely to the averaging done in [13j , 
Equivalently, one could work in Fourier transform space and multiply the Fourier transfori 

by W[j{k), where Wi,{k) interpolates between 1 (for k < 1/L) and for k ^> 1/L. 
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derivatives of short wavelength parts can pick up the factors of 1/e described above, but 
derivatives do not increase the size of long wavelength parts. 

Our framework can be straightforwardly applied to cosmological Newtonian gravity. It 
is natural in this case also to impose the additional conditions that velocities are suitably 
"small" and time derivatives of quantities are correspondingly small compared with space 
derivatives at small scales. Specifically, the sizes we assign^] to the short wavelength part of 
the quantities (ipisr, v l N , 5n) of the previous section are given in Table [fl On the other hand, 



Quantity 


Order 


Wn 


e 


V N 


e l/2 


AS) 
N 


er 1 


do 

d l 


e- 1 



TABLE I. Small scale order counting for Newtonian quantities 



the long wavelength part of all of these quantities and their space and time derivatives 

are assumed to be 0(e). It should be noted that certain products of short wavelength 

quantities can have 0(e) large scale average. In particular, long wavelength averages of 

nonlinear quantities corresponding Newtonian potential energy^, kinetic energy, and linear 

momentum enter the perturbation equation for 7^ . 

Our aim here is simply to use the framework of [1] as a "counting scheme" in powers 

of e to see how close our dictionary comes to producing a solution to Einstein's equation. 

Specifically, we assume that we have been provided with a Newtonian cosmological solution 

where the "sizes" of quantities correspond to Table [fl To complete our "counting scheme," 

(s) 

we must also assign an e-order to B\ . Since Bi is obtained by solving the Poisson equation 
(I2.44p and, according to Table fl], the source term is of order e~ 1//2 , we assign B\ s ^ the order 
e 3/2 . 

Having assigned e-orders to all quantities, we may ask the following question: If we 
substitute the Newtonian solution into our dictionary (I2.40p - fl2.44p to produce a spacetime 
11 The orders we assign to the quantities in Table U correspond to the post-Newtonian orderings of Futamase 

n 

and Schutz up to a rescaling of the spatial coordinates, but they differ from the post-Newtonian 



orderings of Oliynyk 

12 A priori, (5n^n) = 0(1), but due to the fact that 5^ is bounded below by —1, in fact (Sni^n) = O(e); 



see the lemma of section 



II of fl. 
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metric g a b and dust stress-energy tensor T a b, how close does (g a b,T a b) come to satisfying 
Einstein's equation? 

If we are to have confidence that the dictionary is producing a good approximation to a 
solution to Einstein's equation, we would want Einstein's equation to be solved to at least 
0(1) in e. This is a nontrivial requirement, since there are individual terms, such as 5, in 
Einstein's equation that are 0(1/ e) in our counting scheme. As we shall see below, the 
dictionary (I2.40p - fl2.44p solves Einstein's equation to 0(1/ e) but fails to yield a solution 
to Einstein's equation at 0(1) in e. Nevertheless, we will then show that we can make 
further small corrections to the metric so that Einstein's equation does hold to 0(1). As 
we shall see, these metric corrections should be 0(e 2 ) at small scales and therefore should 
be negligible. If so, our original dictionary fl2.40p - fl2.44p should be producing an accurate 
relativistic cosmology in terms of its description of the metric and matter distribution on 
small scales. 

In addition, if the dictionary is to be trusted for its description of large scale structure- 
including on scales comparable to (or larger than) the Hubble radius — we would want Ein- 
stein's equation to hold to at least 0(e) at large scales. As we shall see below, even after 
we have made the necessary corrections to the metric so that Einstein's equation is satisfied 
to 0(1) at small scales, Einstein's equation will fail to hold to 0(e) at large scales in our 
counting scheme. We will therefore make further large scale corrections to the dictionary 
so that Einstein's equation holds to 0(e) at large scales. As we shall see, although these 
corrections are formally of order e, they would be expected to make negligible corrections 
to ordinary linearized perturbation theory at long wavelengths. If so, our original dictio- 
nary ( I2.40p - fl2.44p should be producing an accurate relativistic cosmology in terms of its 
description of the metric and matter distribution on large scales. 

The above corrections provide us with an improved dictionary that incorporates the dom- 
inant general relativistic corrections to fl2.40p - fl2.44p . Although the improved dictionary is 
undoubtedly far more precise than would be needed for most applications, it is important 
as a matter of principle to know that corrections can be made so that Einstein's equa- 
tion holds to 0(1) at all scales and at 0(e) on large scales. Furthermore, for any given 
Newtonian cosmology, the correction terms appearing in the improved dictionary can be 
calculated straightforwardly, and their size should give a reliable indication of the accuracy 
of the original dictionary fl2.40p - fl2.44p . If, as indicated above, these correction terms are 
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negligibly small, then the Newtonian cosmology should provide — via the original dictio- 
nary flZHBD - flgHgp and/or its abridgment fl2T46D - fl2T48l) or simplification f l2T49|) -f l2~5T]) — an 
excellent description of what is predicted by general relativity. 



A. Solving Einstein's equation to O(l) 

1. How well are Einstein's equation solved by the original dictionary? 

Appendix |A] presents the calculation of Einstein's equation for the metric ( I2.36P and 
stress-energy tensor (I2.37l) - (l2.39p . keeping all terms that could potentially contribute to 0(1) 
as well as all terms that could potentially contribute to 0(e) at large scales. Given a New- 
tonian cosmological solution (^n,v 1 n ,5n), we substitute it into the dictionary fj2.40p — (I2.44j) . 
and substitute the result into Einstein's equation, freely using the Newtonian equations to 
simplify the resulting expressions. Equation ( 1A2j) yields 

3 r 2 

V 1 1 3 ' 3 



— <j --d'diipN - -ip N d l diip N - d l i) N dii) N \ - 8np Q {-5 N - (1 + S N )v % N v m } + o(l) 



= ^ ^~ipNd l diip N - d l ip N diip N ^ + 87rpo(l + S N )v l N v Ni + o(l) , (3.3) 

where we used the Poisson equation f!2.23j) for tp^ in the second equality. Since the quantities 
ipNd i di r ipN, d l ipNdiipN and Snv 1 n vni are each 0(1) in our counting scheme (and these terms 
do not cancel), we see that this component of Einstein's equation is not satisfied to 0(1). 
Equation ( 1A3j) yields 



G ^) - 8ttT° = \ {-dii> N - -dtf N \ - 87rpo [(1 + S N )v m ]\ s + o(l) . (3.4) 

However, from the Newtonian equations of motion f!2.23|) and (I2.24p . it follows that 

d'd^N + -tfdtfN = -A7T Po a 2 d l Ul + 6 N )vU . (3.5) 
a 

On the torus this may be integrated, giving 

diipN + -diip N = -4irp a 2 [(1 + 5 N )v Ni ]\ s . (3.6) 

(X 

Thus, we obtain 

G\{g)-^T% = o(l), (3.7) 
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i.e., these components of Einstein's equation are satisfied^ to 0(1). 



Finally, from eq. ( ]A4j) we obtain the space-space components, 

= \ {2$ N - 4tfj N d k d k ij N - 3d k tfj N d k ^ N } 5'j + \ {4^ N d l d^ N + 2d^ N d^ N } 

CL v J CL 

+ i { d ^ j + dj ^) ~ 87rpo(1 + 5n ^ nVn ^ + o(1) • (3 - 8) 

Thus, these components of Einstein's equation are not satisfied to 0(1). 

2. Corrections to the dictionary needed to solve Einstein's equation to O(l) 

We will now show that all components of Einstein's equation can be satisfied to 0(1) by 
making the additional corrections Xi £> an d Jij to the spacetime metric as follows: 

A = ^ N + X + £, (3.9) 
H L = -ip N - X , (3.10) 
hij=jij, (3.11) 

with x^ S \ an d Jij au 0(e 2 ). However, we do not make any modifications to the original 
dictionary expressions for v\ 5, and B\ i.e., we continue to use 



(i + ^y = (i + ^^ + s*)- (i + 5 N y N , (3.i2) 



5 = 5 N - 



. N 2 

a\ a ■ 

- } 1f>N + ~1pN 

a a 



4w 2 i „ ; - , , (3-13) 

d j djB l = -16 W 2 ((1 + 8 N )v% - (1 + 5 N )v^j | . (3.14) 

In particular, it should be emphasized that no additional corrections are made to the matter 
distribution variables 5 and v l . 

We have already seen that the original dictionary solved the time-space components of 
Einstein's equation to 0(1) and it is not difficult to see that these equations continue to hold 
with the above revisions. Thus, to solve Einstein's equation to 0(1), we need only consider 
the space-space components flA4j) and the time-time component (IA2I) . To solve (1A4|) . we 



13 In fact, the precise forms of (|2.4ip and ()2.44|l were chosen so that no further corrections to v l and B l 
would be needed to satisfy Einstein's equation to O(l). 
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note that we can uniquely decompose any symmetric tensor field on a 3-torus with flat 
metric 5ij and flat derivative operator <9, as 

Etj = U6 tJ + d^V - U k d k V + 2d {i W j} + Xij , (3.15) 

with d l Wi = 0, d t X i j = 0, and X l i = 0. This defines the decomposition of E^ into its scalar 
(U,V), vector (Wi), and tensor (Xy) parts. Thus, we can solve an equation of the form 
E^ = by separately solving its scalar, vector, and tensor parts. To begin, we take the 
double divergence of the traceless part of ( 1A4I) . We obtairQ 



3a 2 ' J ^ ' a 2 
1 d l d l {4ij N d j d j ij N + 2d J ij N d j ^j N } 



3a 

= totpodid* [(1 + S N y N v Nj ] - ^-p d% [(1 + 5 N )v k N v Nk ] + o (Jjj . (3.16) 

Here we have dropped terms which are o(l/e 2 ), since we have taken two spatial derivatives 
of an equation that we wish to satisfy to 0(1). We can solve (I3.16P to the desired order by 
defining £ to be the solution to the following double Poisson equation: 

d%d j dj^ = 3did j {2^ N d l d^ N + d l ^ N dj^ N } - {2^ N d j dj^ N + d j ^ N d^ N } 

-127ip a 2 did j [(1 + 6 N )v l N v Nj ] + 4irp a 2 d% [(1 + S N )v k N v Nk ] . (3.17) 

A solution for £ exists on a torus because the source term is a divergence and therefore 
has no spatially constant piece. This solution is unique up to a spatially constant function 
of time, which we fix by requiring that its spatial average, £, vanishes. Since the double 
divergence of the traceless part of (1A4I) has now been solved to 0(l/e 2 ), the scalar part of 
the traceless part of (1A4I) should now be solved to 0(1), as desired. Note that since the four 
spatial derivatives applied to £ yields a quantity that is 0(l/e 2 ), the short wavelength part, 
of £ should be 0(e 2 ), so our assumption that £^ is 0(e 2 ) is self-consistent. 
Next, we show that, with this choice of £, the trace of ( 1A4I) also is satisfied to 0(1). 
Substituting the revised dictionary (13.9p - (13.14j) into the trace of (1A4I) . we find that we must 
satisfy 



a 2 Y 3 *"* T " J a? 

87rp (l + S N )vifV m + o(l) . (3.18) 
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In substituting for we have neglected some terms proportional to 1/(1 + 5m), which can in fact be 



quite large in low density regions. When one makes a uniform momentum correction — (1 + 5j^)v l N 
this corresponds, in a low density region, to a very large velocity correction which is unphysical. If such a 
situation were to occur, then a fix would be to tra^2fcr some of this momentum to a higher density region. 



To see if this equation holds, we take its Laplacian. The double Laplacian of £ will then 
appear, and we can substitute for this quantity using f)3.17p . Since we want to solve ( 13.18)) 
to 0(1), and each spatial derivative increases the small scale order by a factor of 1/e, we 
wish to solve the Laplacian of (13.181) to 0(1/ e 2 ), so the equation we wish to solve is 

^ [2d%4> N - 2d 1 {drtN&d&x) } = 24np 8^d t [(1 + 6 N )v%v Nj ] + o (±\ . (3.19) 

However, using the Newtonian equations (I2.23p - fl2.25p . as well as the Friedmann equations 
for the Newtonian background, one can show that 



d l di <J 2$ N + 6-tp N + 
a 



a J \a 



4>N 

= 8 W 2 <9^ [(1 + 5 N y N v j N ] + 2di (d j d^ N d^ N ) , (3.20) 

so f )3.19p is indeed solved to the desired order. We have thus fully solved the scalar parts of 
the space-space part of Einstein's equation to 0(1). 

Next, we consider the vector part of (IA4[) by taking its divergence, using the fact that 
the scalar parts have already been solved. We obtain 

^8 3 & + i# {Afa&dtfN + 2d^ N d^ N } I 

= 8np 8 j [(1 + 5 N )v%v Nj \ \ v + o Qj . (3.21) 

Using the Newtonian Euler and mass conservation equations, along with the definition of 
Bi, one can show that 

^ [cPfyB* + 2°-&d j B^ + ^ {A^ N d l d^ N + 2d^ N d^ N )l 

= tocpotf [(1 + 5 N y N v Nj ] I , (3.22) 

so equality does hold for the terms explicitly written in f l3.2ip . Thus, the vector part of 
(fAil) is satisfied to 0(1). 

The tensor part of (1A4|) is all that remains of this equation. We obtain 

- \d k d k f 3 + \ {tyN&djfa + 2d^ N d^ N }\ t = 87r Po (l + 5 N )vifV Nj \ t + o(l) , (3.23) 

where \ t denotes the tensor part of a quantity in its decomposition f )3.15p . To solve this 
equation to the desired order, we define jij to be the solution of 

8 k 8 k f j = {87rp a 2 (l + S N y N v Nj + A^N&dj^N + 28 i ^ N 8 j ^ N 

-87rp o 2 (l + M<^ + 28^x8^} ^ , (3.24) 
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where the overline denotes spatial average (see (I2.45P ). The terms with the overline in (I3.24p 



, as is necessary in order to be 
0(e), we will satisfy (I3.23P to 



have been added in so that the source has vanishing integra 
able to solve the Poisson equation. Since these terms ar 
the desired order by choosing to solve f)3.24p . We fix the ambiguity in by requiring 
Jij = 0. Since two spatial derivatives applied to jij yields a quantity that is 0(1), the 
short wavelength part, of jij should be 0(e 2 ), so our assumption that is 0(e 2 ) is 
self-consistent. We have now solved ( 1A4|) to 0(1). 

Finally, we consider the time-time component of Einstein's equation. Substitution of the 
dictionary into (IA2[) yields 

3 f 2 2 8 

— <^ --d'diipN - -d'diX - -^Nd l diifj N - &il) N diil) N 
a z { 3 3 3 

= 87rpo [S N - (1 + 6 N y N v Nl ] + o(l) . (3.25) 

Using the Newtonian field equation (12.231) . we obtain 
3 f 2 



i d l diX - -^ N d l diij N - d'ipNdiipN } = -8np (l + 5 N )v l N v Ni + o(l) . (3.26) 
a z { 3 3 J 

We define x t° be the solution to the Poisson equation 

3 

d l diX = -Aip N d l d^ N - -d l ip N diip N + 4vra 2 po(l + 5 N )v l N v Ni 



--ffipNdiipN - Ana p (l + d^v^vm , (3.27) 

with x — 0. Since two spatial derivatives applied to x yields a quantity that is 0(1), the 
short wavelength part, x^ S \ °f X should be 0(e 2 ). 

Thus, we have shown that Einstein's equation can be solved to 0(1) by making the 
corrections fl3.9p - fl3.lip to the original dictionary, where x, £, and jij are given, respectively, 
by (I3.27p . f)3.16p . and f)3.24p . Although it is extremely important as a matter of principle 
that such corrections can be made so as to obtain a solution to 0(1), we expect that these 
corrections will be negligibly small compared with ipjy. 

B. Improving the solution to 0(e) at large scales 

In the previous subsection we obtained a solution to 0(1). However, as previously stated, 
if our dictionary is to be trusted for its description of large scale structure — including on 
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See footnote [T2l 
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scales comparable to (or larger than) the Hubble radius — we want Einstein's equation to 
hold to at least 0(e) at large scales. Within the context of ordinary perturbation theory, 
this corresponds to solving the linearized perturbation equation. Within our generalized 
perturbative framework, this corresponds to solving the generalized linearized perturbation 
equation (87) of The difference between these, as noted earlier, is that long wavelength 
averages of products of small scale quantities enter into the generalized linearized equation. 

It is easy to check that, even with the corrections fl3.9p -f l3TTTj) . our dictionary does not 
produce a solution to 0(e) at long wavelengths. Therefore, we will need to make the following 
additional long wavelength corrections to our metric and matter variables: 



.4 


= ii N + x + Z + x +~, 






(3.28) 


H L 


= -4>n - x - X , 






(3.29) 


1 + 5 N )vi 


= (1 + 6 N )(v m + Bi) - (1 + 5 N )v Ni 


+ P i} 

V 


(3.30) 


5 


-s N 3 


f a\ 2 , a ■ 

- Vn + ~1f>N 

\a J a 


+ A, 


(3.31) 












(3.32) 



No additional long wavelength correction is needed for B^. Here, the quantities 5, X, Pj, A, 
and are assumed to be 0(e) and to have vanishing short wavelength part. Hence, they 
do not contribute to Einstein's equation to 0(1) and thus do not spoil the solution obtained 
in the previous subsection. 

Our strategy is to apply the averaging operator (•) (see eq. (13.11) ) to Einstein's equation, 
and to choose the above new correction terms in order to obtain a solution to 0(e). For 
our the calculations below, it is useful to note that the averaging operator (•) commutes 
with differentiation. Note also that since ipN = 0(e), we clearly have (ip%) = 0(e 2 ), and, 
consequently, we have 

(d^ N d^ N ) + & N d%i; N ) = \d%{^ 2 N ) = 0(e 2 ) . (3.33) 

Thus, we may freely "integrate by parts" to set (ipMd % diip^) = —(d l ipNdiipN) in our calcu- 
lations. 

As before, we begin with the double divergence of the trace-free part of the space-space 
components (1A4j) of Einstein's equation. Substituting our new dictionary fl3.28p - fl3.32p . ap- 
plying the averaging operator (•), and using the equation obtained by applying the averaging 
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operator to (I3.17P to simplify the resulting expression, we obtain 

- ^d%&d 3 ~ = o(e) . 
Thus 5 can only have a spatially constant part, i.e., 



(3.34) 



(3.35) 



where S may be an arbitrary function of r. Examining the scalar homogeneous parts of the 
metric, 

ds 1 = a 2 (r) [-(1 + 2X + 2E)dr 2 + (1 - 2X)5 ij dx i dx j ] , (3.36) 
we see that S corresponds to gauge freedom in the choice of time coordinate. We fix this 



freedom by setting 



-2X 



(3.37) 



corresponding to using conformal time. 

Next, we consider the trace of (1A4j) . Substituting from the dictionary, applying (•), and 
using H = — 2X, we obtain 



3 r 2 



a 2 3 



<m<£) + 2 (4> N ) + (x)+X)+2- 3(^) + 3(x) + 3X + (0 - 2X 



+ 



49 T - +2 - 



= 87rp ((l + M^Wi) + o(e). (3.38) 

As before, we take the Laplacian of this equation and substitute the average of (13.171) . 
obtaining 

- 2 &d 3 (2 ($ N ) + (x) + x)+2- (3(ip N ) + 3(x) + 3X + (£) 



+ 



4<9 T - +2 - 



(& N ) + (x) + x + (0) > - ^(3V;v^iv> 



247rp ^^'((l + (Jjv) «Wi> + o(e) . 



(3.39) 



Simplifying further using the average of f)3.20p and then inverting the Laplacian, we obtain 

2" 



a • a— 
2X + 6-X-4-X + 

a a 



4a. 



+ 2 



a a 



45L 



+ 2 



(X-X) 

«x) + (0) 



1 



d^NdiipN 



mp a 



1 + S N )v % N v Ni + o(e) . 



(3.40) 



26 



Here, the constant of integration was determined by requiring consistency with (13.381) . Thus, 
the trace of (1A4|) is satisfied to 0(e) at long wavelengths provided that X satisfies this second 
order ordinary differential equation in time. The "scalar parts" of the long wavelength part 
of (IA4|) have now been satisfied to 0(e) 

Using (13.221) . it is not difficult to see that the long wavelength part of the divergence 
of ( 1A4I) is solved to 0(e) at large scales, without any need for further corrections. Thus, 
the "vector part" of (IA4|) has been satisfied to the desired order at long wavelengths. Only 
the "tensor part" of (1A4j) remains. Substituting the dictionary and applying the averaging 
operator (•), we obtain 

\\j i j + -j i j -d k d k J i . 

= ~\ {(./',) + ~(fj) - d k d k (fA + - 2 (d^ N d^ N )\ t + 8vrp ((l + 6 N )v%v Nj )\ t + o(e) 

(X I Gi I d 



= {OV + + 4^^ +8irp {l + 5 N )v i N v Nj +o(e), (3.41) 

where we used the average of (13.241) in the second line. Thus the tensor part of flA4j) is 
solved to 0(e) at large scales provided that Jij solves this wave equation. This completes 
the solution of the long wavelength part of flA4j) to 0(e). 

Next, we consider ( 1A3I) . Substituting from the dictionary, applying (•), and taking the 
divergence, we obtain 

-- 2 d% l$ N ) + x + ( x ) + - (& N ) + x + ( x ) + (0) 

= torpo& ((1 + 5 N )v m + Pi) + o(e) . (3.42) 
Using the average of (13.51) to simplify this expression, we obtain 

- \x + ( X ) + - (X + (x) + (0) ) = torpdtPi + o(e) . (3.43) 



We solve this equation by setting P l to be 



' r d i (x + (x) + -(X+( X ) + (0)) ■ (3-44) 



ATrp a 2 \ ' a 

This satisfies the "scalar part" of (1A3I) to the desired order at long wavelengths. It is easy 
to check that the vector part of (1A3I) is also satisfied without the need for any further 
corrections. 



27 



Finally, we consider the remaining component of Einstein's equation, the time-time com- 
ponent (1A2|) . Substituting and averaging, we find that this equation is satisfied to the 
required order by making the density correction 

iira z p 67rp a z { 3 



A = -((1 + 6 N y N v Nl ) - -——(^diTpff) - - — 2 { --d l d t (X + (x» 



+2- a (X + <*)) + 2 Q 2 (X - 2X + ( X ) + (0) 



5 3 r 2 

-(1 + 6 N )v* N v m - - — —d^NdiipN - <^ --d l diX 

S7ra z po 57rpoo I 3 



+2^ (X + ( X )) + 2 Q (X - 2X + (x) + (0) I • (3-45) 

Here, the average of (13.271) was used to get the second line. 

Einstein's equation has now been fully solved to 0(1) everywhere, and to 0(e) at large 
scales. All of the quantities appearing in our dictionary are uniquely determined by the 
Newtonian solution, except for X and J^-, which obey second order differential equations in 
time. The degrees of freedom associated with X correspond to the long wavelength degrees 
of freedom present in the dust matter sector in ordinary linearized perturbation theory. It 
would be natural to fix X by requiring that A and Pi vanish at an initial time 16 !. The degrees 
of freedom associated with Jy correspond to the presence of long wavelength gravitational 
radiation. 

Finally, we consider the magnitude of the additional long wavelength quantities 5 = — 2X, 
X, Pi, A, and that we have just obtained. The equations for these quantities involve 
terms of the form {diifj^djipN) and ((1 + 8j^)v m v Nj) as we U as (£)> (x)> an d (jij), which 
themselves are sourced by terms of the form (diipNdjipN) and ((1 + 6jf)vmVNj)- Thus, the 
additional long wavelength quantities appearing in our new dictionary (I3.28p - fl3.32p should 
have a magnitude of the order of the Newtonian potential energy and kinetic energy of 
the dust matter. Although very small, the homogeneous (i.e., spatially constant) part of 
these quantities provides the dominant correction to the background FLRW dust cosmology. 
We compute these corrections explicitly in Appendix [Bl However, the long wavelength cor- 
rections at finite wavelength are sourced by large scale inhomogeneities in (diifj^djipN) and 
((l + 5^)vmVNj)- These terms should be extremely small as compared with, say, {Sn)- Thus, 
16 Since Pi vanishes identically, this does not fix the spatially homogeneous part, X, of X. An additional 
condition on X will be imposed in Appendix [B] . 
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the long wavelength corrections we have obtained in this subsection should make entirely 
negligible contributions to Newtonian large scale structure. 



IV. SUMMARY 



Combining all of the results of the previous section, we have the following Oxford dictio- 
nary for translating a Newtonian cosmological solution (ipN, v l N , 6~n) to a general relativistic 
spacetime metric f!2.36j) and dust stress-energy (12.371) : 



A = ^ N + X + Z + X -2X 
H L = -ip N - x ~ X , 



1 + 5 N )vi = (1 + 5 N )(v m + B,i) - (1 + 5 N )v N 

2 

.-> / II \ 

5 = 5 N - 



4iip a 2 

h%j Jij ~\~ J%j ■ 



a \ , a ; 

- ] IpN + ~WN 

a a 



+ Pi, 
+ A, 



(4.1) 
(4.2) 
(4.3) 

(4.4) 
(4.5) 



Here, the quantities Bi, £, x, j^, P\ and A are given by 



&>djB i 



pi 

A 



167rp a" (1 + 8 N y N - (1 + tf^v 



iV 



d^^dji = 3d t d J ^ti&drft, + tf^tidrfx} - {2^ N d j d j ^ N + &>^ N d^ N ) 

-127Tp a 2 did J [(1 + 5 N )v l N v Nj ] + ATtp^d'di [(1 + 5 N )v k N v Nk ] , 
3 

d l diX = -4-0 n9' diip N - -d'tppfdiipN + 4?ra 2 po(l + S N )v l N v m 



-^ipNdiipN - Atto p (l + S N )v l N v Ni , 
d k dkfj = {87rp a 2 (l + S N y N v Nj + Ai) N d l dy N + 2d i ip N d j ip N 



-8np a (1 + 8 N )v l N v Nj + 2d i i> N d j 'ijj N 



1 d i (x + ( x ) + -(x + (x) + (0) 



4:7ipoa 



-(1 + 8 N )v z N v Ni - 



]ira 2 po 



d^NdiipN 



+2-[X+( X })+2 



mpoCL 

(X-2X+( X ) + (0) 



--d'diX 
3 



(4.6) 
(4.7) 

(4.8) 

(4.9) 
(4.10) 



(4.11) 
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with Bi = £ = x — Jij ; — 0. The quantities and X satisfy the differential equations 



and 



J% + -j\ - d k d k J\ 



(fj) + -(fj) \ + — d^Nd^N + 87rpo(l + S N )v l N v Nj 



a ■ a— 
IX + 6-X - A-X + 

a a 



id, I 5 +2(5 

a / \ a 



-2(x)-6-(x)-2-(0- 
a a 



idA a - +2 (i 



((x> + (0) 



i- 



irpoa 



--O'c.xili/x + o 



(f + 5 N )v l N v Ni 



(4.12) 



(4.13) 



Like the Oxford English Dictionary, the above dictionary should be far more detailed and 
precise than needed for everyday use. Nevertheless, it may be comforting to have it on 
one's bookshelf in case the need does arise. Furthermore, as a matter of principle, it is of 
importance to know that a dictionary of this accuracy — namely, solving Einstein's equation 
to 0(1) on all scales and to 0(e) on large scales — can be constructed without running into 
inconsistencies. 

Our main purpose in obtaining the complete dictionary ( |4.ip - (j4.5p was to evaluate the 
accuracy of the original dictionary (I2.40p - fl2.44p (as well as its abridgement (I2.46p ~ (j2.48p and 
simplification (I2.49p — (I2.5ip ). We have argued that for a Newtonian cosmology that satisfies 
ipN 1 and | i^l 1 but may have 5n ^> 1 at small scales, all of the additional terms 
appearing in ( 14. Ill — ( 1431) as compared with ( I2.40p -f f2.44p should be negligibly small. Whether 
or not this is actually the case for any given Newtonian cosmology can be determined by 
computing the quantities £, x, Jij, P\ A, X, and given by (I4.7p - (l4.13p . If these quantities 
are indeed negligibly small, then one can have confidence that the Newtonian cosmology is 
accurately representing a general relativistic spacetime via the original dictionary ( 12.4011 — 
(I2.44p . If, in addition, B l is negligibly small (see (I4.6P ). then one is similarly justified in 
using the abridged dictionary ( 12.4611 — (12.4811 . These statements remain valid even in cases 
where the Newtonian cosmology is describing phenomena on scales comparable to or larger 
than the Hubble radius. 
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Appendix A: Einstein's equation 

For the metric ansatz ( 12.361) . we write down the various components of the perturbed 
Einstein equation, 

G*M - G^(9 m ) = Svr (n - ^ (0)M „) , (Al) 
keeping all terms which are 0(1) at small scales or 0(e) at large scales in our counting 
scheme. We introduce the notation o(l; e) to denote a quantity which is o(l) at small scales 
and o(e) at large scales. The li v = ° equation reads 
3 I „d • 2 „ „ /a x 2 



-2-F L + -d'c^ + 2 - A - -H L d l d t H L - d l E L d l R L 
a z \ a 3 \ a / 3 

= Svrpo {-5 - (1 + Sy( Vl - Bi)} + o(l; e) , (A2) 
the ^ v = ° i equation is 

^ \d t H L - %A\ - -^&d s Bi = 87rpo(l + 6)(Vi - B { ) + o(l; e) , (A3) 

and the M „ = * • equation is 

' W(# L + A) -2tf L - 4^ + 2^4 + 4^ + 2 ("Y^M' 



a 2 I a a \a / \a 

+^ {-AH L d k d k H L - 2Ad k d k A - 2H L d k d k A - 2d k H L d k H L - d k Ad k A] 5) 

—z&dAH L + A) + 4 UE L d i d 1 H L + 2H L d i 8 i A + 2A#&j4 + <9Mcl+L 
+d i Ad j H L + &H L djA + Z&H L djH L } 

+-^ (9*4 + 6^ + 2-d%- + 2-g^I + -4 U\ + - ^SlfcA*- j 
2a 2 t a a J a 2 [ 7 a J 7 J 

87rpo(l + S)v\vj - Bj) + o(l; e) . (A4) 



Appendix B: Modified background metric 

In this Appendix we compute the homogeneous part of the metric and matter distribution 
as given by our final dictionary (I4.ip - fl4.5p . These can be viewed as providing the dominant 
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corrections to the background cosmology produced by small scale inhomogeneities. The 
relevant equations can be obtained by taking spatial integrals of the equations of section HVl 
We find that the spatially homogeneous parts of the metric components are given by 



A 
Hl 



-X 
-X 



hij J ij ; , 



(Bl) 
(B2) 

(B3) 



as well as B l = (see footnote [5]) . Thus, the homogeneous part of the metric takes the form 



ds 2 = a 2 (r) [-(1 - 2X)dr 2 + ((1 - 2X)5 ij + J^d^dbP] . 



(B4) 



We also have 



5 = A 



and 



(1 + S)(vi-Bi) = 0. 
In addition, the quantities X, A, and J ^ satisfy 



2X + 2-X 

a 



a } \a 



X 



-d^jsfdi^N + ^a 2 po(l + 5 N )v l N v Ni , 



(B5) 
(B6) 

(B7) 



fB8) 



A = -(1 + S N )v l N v Ni 



^7ra 2 po 



2 a -X-2( a -\ X 



and 



1 Jy a^f 



T d<"i(; N djip N + 87rp (l + S N )v l N v Nj 



(B9) 



(BIO) 



It is clear that the metric perturbation given by X can be interpreted as taking one to a 
new FLRW spacetime, with scale factor 



a( T ) = a(r)(l-X) . 



(Bll) 



We now derive modified Friedmann equations for a. To linear order in barred quantities, we 
have 



Ida 1 da 
adr adr 
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X 



(B12) 



so 



1 dd\ 2 /l<ia\ 2 2®j( 
adrj \adrj a 



" l>[)(r ( \ + A + (1 + 5 N y N v N ^ + ^ViV^iV + ■ (B13) 



3 V v JV "V 3 3 

Here, we made use of the Friedmann equation for a as well as eq. ( 1B9I) . Similarly, we have 

<i / 1 da \ d ( 1 da \ +> 
dr \ddr J dr \adr ) 



3 



(l + A + 2(1 + 5ivK^) - ^ N d^ N + -|- . (B14) 



To put these equations in a more recognizable form, we note that for dust matter, 
V a (/7u a ) = 0, so the integrated flux of pu a over a Cauchy surface £ 

M = - J pu d e dabc (B15) 

is a constant, i.e., independent of E. M is often referred to as the "total number of baryons"; 
in an iV-body simulation, it would correspond to the total number of particles in the simu- 
lation. Evaluating the right side of (IB 151) . we obtain 

M = Po a 3 ( 1 + A + 1(1 + 5 N y N v m + -J^d^ N d^ N - 3x] . (B16) 



It is natural to use our freedom in choosing initial conditions for a solution to (JB8J) to require 
M = Ao = poa 3 , so that the total number of particles is the same as in the background 
spacetime. This condition yields 

-r- 1- 



A = --(1 + 5 N y N v Nl - ^-—^d^ N d^ N + 3X . (B17) 



2 47rp a 
Combining this equation with (1B9I) . we obtain 



1^ =— ^ : 1 -T— : — ; — — 3 I d^-r f (1^ ^ 



= --(1 + 5 N )v%v m + -—-fr^d^x -3X- -— =- i2-X-2[-) X } , (B18) 
2 8na z p 8na 2 p | a \a/ 1 

whose time derivative is flB8j) . 

We define the average particle number density p relative to a by 

pa 3 = Af = Af = p a 3 . (B19) 
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In terms of p, the Friedmann equations become 

lda\ 2 8npd 2 / 1-. ^--^ 1 — — — — — \ Ad 2 /T ^„ N 

-aTr) = -f" I 1 + 2* 1 + " g^P**"**" ) + T ' <B20) 

" ^ ^ ' 1 + + WW- - A^MA + ^ . (B21) 



From these equations, one can read off the effective energy density and pressure, including 
the contributions from small scale inhomogeneities, 



Pes = p ( 1 + -(1 + S^v^VNi - Q _^ 2 d i ^ N d i ilj N ) , (B22) 



U 1 



Pes = P o(l + S N )vifV N i - ^y ^jfdi^N . (B23) 
\o I'm pa J 

The correction terms in p e fj correspond precisely to averaged gravitational potential energy 
and kinetic energy, as expected |l(J, ll3j. They can be interpreted as renormalizing the 
proper mass density p to an "ADM mass density" p e s- For virialized systems, the correction 
terms in P e fj cancel, as pointed out in fl3 |. Thus, we see that the corrections resulting from 
X and A correspond to modifying the FLRW background to a new FLRW spacetime with 
small corrections to the average effective mass density and pressure that arise from small 
scale Newtonian gravitational potential energy and stresses as well as small scale kinetic 
motions. 

The remaining corrections due to Jj 3 - perturb one to an anisotropically expanding Bianchi 
model. It can be seen from (IBlOj) that anisotropics in the spatial average of the Newtonian 
stresses and/or kinetic motions must necessarily induce an anisotropic expansion of the 
universe. However, we would expect these effects to be extremely small. 
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